The Role of Chaos in One-Dimensional Heat Conductivity 
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We investigate the heat conduction in a quasi 1-D gas model with various degree of chaos. Our 
calculations indicate that the heat conductivity k is independent of system size when the chaos of the 
channel is strong enough. The different diffusion behaviors for the cases of chaotic and non-chaotic 
channels are also studied. The numerical results of divergent exponent a of heat conduction and 
diffusion exponent (3 are in consistent with the formula a = 2 — 2/13. We explore the temperature 
profiles numerically and analytically, which show that the temperature jump is primarily attributed 
to superdiffusion for both non-chaotic and chaotic cases, and for the latter case of superdiffusion 
the finite-size affects the value of /3 remarkably. 

PACS numbers: 44.10.-|-i, 05.45.-a, 05.70.Ln, 66.70.+f 



I. INTRODUCTION 

The low dimensional microscopic dynamics of heat 
conduction has been an attractive question since early 
year last century. Much more attention has been payed 
to this problem in the last two decades due to the 
dramatic achievement in the application of miniatur- 
ized devices Q, IS IE ^ IE IS which can be described 
by ID or 2D models. More and more numerical cal- 
culations are focused on the minimal requirements for 
a dynamical model where Fourier's law holds or not 
HlilSIIillililllSlilllllll- a convergent heat 
conductivity was shown in ding-a-ling model lajp w hich 
is chaotic. The studies on Lorentz gas model [iflllll (the 
circular scatters are periodically placed in the channel) 
of which the Lyapunov exponent is nonzero gave a finite 
heat conductivity which fulfils the Fourier law explicitly. 
Hence, chaos was ever regarded as an indispensable fac- 
tor to normal heat conduction. Whereas, the FPU model 
[l^ IT^ indicated that the chaotic behavior is not suffi- 
cient to arrive at normal heat conduction. Recently, a 
series of billiard gas models ^10, J4, .15,, were devoted 
to explore the normal heat conduction of quasi ID chan- 
nels with zero Lyapunov exponent. However, the role of 
chaos in heat conduction has not been well understood. 
Additionally, the exponential stability and instability fre- 
quently coexist in the scatters of real system. Thus the 
model with various degree of chaos deserves further in- 
vestigation from the microscopic point of view, and it will 
be also interesting to explore non-equilibrium stationary 
states and to determine the steady temperature field. 

In this paper, we focus on the quasi one-dimensional 
gas model which is closer to the real system. The scat- 
ters in our model are the isosceles right triangle with a 
segment of circle substituting for the right angle. In this 
case the edges of scatters are the combination of line and 
a quarter of circle. Such a channel is of chaos, which 
indicates exponential instability of microscopic dynam- 
ics. Our paper is organized as follows. In section^ we 
introduce the model and investigate the degree of chaos 
for various channels with different arc-radius and channel 



height. In section lHIl we study the heat transport behav- 
ior and the corresponding diffusive behavior by changing 
the radius of top arc and channel height. In section HVI 
we investigate the non-equilibrium stationary state and 
determine the steady temperature field numerically. We 
also analyze the dependence on diffusion exponent fi and 
system size A'' of temperature profile theoretically. In 
section we discuss the relation between our work and 
others and summarize our main conclusions. 



II. THE MODEL 

We consider a billiard gas channel with two parallel 
walls and a series of scatters. The channel consists of N 
replicate cells of length I and height h, and each cell is 
placed with two scatters as shown in Fig. ^ The scat- 
ter 's geometry is an isosceles right triangle of hypotenuse 
a whose vertex angle is replaced by a segment of circle 
with radius R which is tangential to the two sides of the 
triangle. At the two ends of the channel are two heat 
baths with temperature Tl and T^. Noninteracting par- 
ticles coming from these heat baths are scattered by the 
walls and the straight lines as well as the arcs of the 
scatters in the channel. 
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FIG. 1: The channel with N replicate cells. Here, / = 2.2, 
a = 1.2, h changes from 1.0 to 0.27, and R from to 0.848528 
to ensure a quarter of circle always. 

For such a channel, the degree of chaos can be charac- 
terized quahtatively by Poincare surface of section (SOS) 
|r^ . Suppose we take out one unit cell from the channel 
and close the two ends by straight walls. Then the prob- 
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lem becomes a billiard problem. A particle moves within 
the cell and makes elastic collision with the mirror-like 
boundary. We investigate the surfaces of section (s, Vr) 
under different initial conditions, s is the length along 
the billiard boundary from the collision point to the ref- 
erence point. Vr is the tangential component of velocity 
with respect to the boundary at that point. The filling 
behavior of phase space shown in Fig. |21 indicates the 
degree of chaos. In case I, it's non chaotic. The surface 
of section is regular and periodic as shown in Fig. ^a). 
As the radius R is increased from (a) to (e), the motion 
becomes more complex and the map becomes dense with 
points except some regular islands. In Fig. [3f), the 
regular parts disappear which indicates strong chaos. 
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FIG. 2: Poincare surface-of-section of the billiard problem. 
The billiard starts with an incident angle 0.8 and unit velocity, 
(a) 7? = 0, /i = 1.0; (b) R = 0.001, h = 1.0; (c) R = 0.015, 
h = 1.0; (d) R = 0.1, h = 1.0; (e) R = 0.848528, h = 1.0; (f) 
R = 0.848528, h = 0.27. 



III. THE HEAT TRANSPORT AND DIFFUSION 
BEHAVIOR 

To study the heat conduction of the model, the heat 
flux is investigated firstly. In calculating the heat flux, we 
follow Ref. ■ For simplicity, the particles from the two 
heat baths are supposed to have definite velocities \/2TY 



and \/2Tr respectively 0|. We consider one particle 
colliding with a heat bath during a period of simulating 
time. The energy exchange (AE)j at the jth collision 
with the heat bath is defined as 



En 



(1) 



where E^ denotes for the energy of particle taken from 
the bath and Ep for that carried in the channel. For M 
collisions between the particle and the bath wall during 
the simulation time t, the heat fiux is given by 



t 



(2) 



As there is one heat carrier in each cell and the channel 
has TV replicas, there are N particles in the whole chan- 
nel. Summing over the heat flux of N heat carriers, we 
have Jn{N) = NJi{N). Meanwhile, the Fourier's law 
reads 



JNiN) 



dT 
dx 



Tl-Tr 

' m ' 



(3) 



where k refers to the heat conductivity which is deter- 
mined by Eqs. 10 and ©, 



(4) 



We consider various cases by changing the radius R 
of the top arc of the scatters to investigate their effects 
on heat conduction. The heat flux of a single parti- 
cle versus system size shown in Fig. |31 are four typical 
cases. Namely, case I: the ■ studies for i? = 0, ft- = 1.0; 
case I: the o for i? = 0.001, h = 1.0; case IE: the A 
for R = 0.848528, h = 1.0, and case IV: the V for 
R = 0.848528, h = 0.27, respectively. The total cell 
numbers are chosen as iV = 20, 40, 80, 160, 320, 640 
and 1280 respectively. After a sufficient long period of 
simulation time, the heat fiux approaches to a constant 
value. Clearly, the value of heat flux decreases with in- 
creasing R for the same size. Remarkably, there is 20 
times difference of heat flux between case III and case 
IV, which indicates that smaller height suppresses the 
heat flux greatly. Thus, it appears that the value of heat 
flux can be adjusted in this way in designing heat-control 
devices. Furthermore, our calculations show that the 
heat flux dependence on N exhibits faint non-linearity 
although the curve looks linear for all cases except case 
IV in the log-log scale . 

In order to observe the deviation from the line, which 
arises from the finite-size effect, we calculate the ratio of 
heat fiux versus system size for various radius. The data 
for the aforementioned four cases are plotted in the right 
panel of Fig. |21 from which one can see that both the 
increasing of system size and of the arc radius bring the 
ratio an upward tendency to the value 4 which ensures 
the Fourier law. In case I where i? = 0, there is only a 
slight increase for the ratio around 2.4. Whereas, the ra- 
tio rises drastically along with the increasing of systems 
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FIG. 3: The heat flux of a single particle versus system 
siz6(Af = 20, 40, 80, 160, 320, 640 and 1280) with the diver- 
gence exponent of heat conductivity q = 0.721, 0.526, 0.101, 
0.009 for four typical cases respectively (left panel). The ra- 
tio of heat flux Ji (A'^) / Ji {2N) for different system sizes (right 
panel) . 



size even if the radius R is merely 0.001 (the case I). 
When R — 0.848528 (the case H), the scatters become a 
full segment of quarter circle. The ratio also rises dras- 
tically at first and gradually after N > 160 in this case. 
In both cases I and IH, it seemly approaches to distinct 
asymptotic values which are all different from that for 
normal conduction. This implies that smaller degree of 
chaos is insufficient to bring about a normal heat con- 
duction although the increasingly chaotic degree makes 
the divergent exponent of heat conduction smaller. In 
case IV, we maintain the scatters at radius R = 0.848528 
and reduce the height h from 1.0 to 0.27. In this strongly 
chaotic case, the ratio fluctuates around the value 4 (dot- 
ted line) which means that Fourier law is obeyed. 

It is known that the normal heat conduction happens 
when a = 0, which indicates the heat conductivity is 
independent of system size, and the anomalous heat con- 
duction corresponds to the case of a > 0. The heat con- 
ductivity K we calculated can be given by k ~ N"' with 
a > despite the heat-flux ratio has a different increase 
in asymptotic value for all cases (except case IV). 

We calculate a at the range of system sizes N from 20 
to 1280 by averaging over many realizations for various 
radius R at fixed channel height h = 1.0, and the plot of 
the dependence of a on i? is shown in Fig. QJa). One can 
see that a descends from 0.721 through 0.526 to 0.092 
if R increases from to 0.848528 for a fixed height h = 
1.0. Clearly, the a descends rapidly for small radius {e.g. 
R ~ 0.001 in case 1) and slowly for larger ones. This 
illustrates that the appearance of arc on the top of the 
scatter suppresses the divergent exponent a drastically. 
If the channel height h for fixed R = 0.848528 is changed 
from 1.0 to 0.27, the a is found to diminish to 0.009. 
Therefore the k appears to be independent of system size 
and the Fourier law holds in this case. 

Since the characteristic of heat transport is found be- 
ing closely related to the diffusion behavior 0, 0, 
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FIG. 4: (color on line) (a) Conductivity divergence exponent 
a versus circular radius R. The ■, refers to the magnitude 
of a for /i = 1.0, R = 0, 0.001, 0.05, 0.1, 0.2, 0.4, 0.6, and 
0.848528; the A for /i = 0.5 and R = 0.848528; the * for /i = 
0.27 and R = 0.848528 has the value of 0.009. (b) The relation 
between /? and a, where the circle is the numerical result, the 
red line is of a = 2 - 2//3 fl^ and the dashed line is the result 
of Ref. (c) Log-log plot of mean square displacement 

{x{t)^) versus time t. The curves from top to bottom on the 
right correspond to cases I, H, HI and IV respectively. The 
ensemble has 10^ particles starting from the center of the 
channel at time t = where a; = with the unit velocity and 
random direction. 



IT^ I20I I , we investigate the diffusion property for the 
above cases subsequently. For a particle starting at the 
origin at time t = and diffusing along x direction, the 
mean square displacement {{x{t) — a;(0))^) characterizes 
its diffusion behavior. For normal diffusion, the Einstein 
relation of {{x{t) — a:(0))^) = Dt holds, where D is dif- 
fusion coefficient. If the mean square displacement does 
not grow linearly in time, i.e., {{x{t) — a::(0))^) = Dt^ , 
we refer to anomalous diffusion. Recently, the connec- 
tion between anomalous diffusion and corresponding heat 
conduction in ID system was discussed hotly [isL 113. l20| . 
We plot the mean square displacement versus time t in 
Fig. ^c) for the aforementioned four cases. Note that 
10^ particles were put at the center of the channel where 
X — with unit velocity and random direction in the 
simulations. The top solid line and the bottom short- 
dash-dot line are precisely straight in the whole Simula- 
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tion period {t = 10^), which correspond to the case I 
(non-chaotic) and case IV (strong chaotic), respectively. 
We obtain f3 = 1.628 for the case I which corresponds to 
a = 0.721, and (3 = 1.001 to a = 0.009 for the case IV. 
Beyond these two cases do the curves keep asymptotically 
linear at large time t with diffusion exponent (3 between 
the values of above two cases. The best fits of the slope 
give (3 = 1.357 which corresponds to a = 0.526 for case 
I and f3 = 1.050 to a = 0.101 for case IE, respectively. 
The relation between divergent exponent a and diffusion 
exponent /3 fits the relation oi a = 2 — 2/f3 proposed by 
Li and Wang in Ref. 0, as is plotted in Fig. ^h). 
Whereas Denisov et al. presented another connection of 
a with /3 on the basis of the Levy walk model • 2i| . More 
details about the origin of the discrepancy between above 
two relations can be found in Ref. [13 . 
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FIG. 5: (color on line) (a), (b), (c) and (d) The PDFs of 
the flight distance \Sx\ between two consecutive collisions for 
above four cases I, I, H and IV respectively. A'' represents the 
system size. Note that the log-log scale is used in (b) and the 
inset of (d). In the latter case (d), Gaussian distribution (red 
dashed-line) is in comparison to the numerical PDF. (e) The 
typical trajectory with periodicity in case of = , ft = 1.0. 

As different diffusion behaviors are likely related to the 
trajectory characteristics of the particle propagation, we 
investigate the PDF ^(|(Sa;|) of the flight distance \Sx\ in 
x-direction between two consecutive collisions with the 
scatters. After a long time for adequate collisions in the 
channel, the PDF for aforementioned fore cases, shown in 
Fig. [Sfa) to (d) respectively, take on completely different 
forms for different cases. In case I, the discrete values of 
probability indicate that the trajectories are abundant of 
periodicity, which is almost alike for larger system size. 



The maximum value of PDF appears when \Sx\ = 0.447, 
and the typical trajectory is plotted in Fig. ]^e) which 
shows explicitly that the parallel passage makes the pe- 
riodical trajectory possible, and the particles are easier 
to propagate along the channel with fewer collisions. It 
is superdiffusion in this case. In case E, only smaller 
system size has the explicit periodicity. With the system 
size growing, the periodicity is destroyed by the collisions 
with the segment of circle time and time. The PDF gets 
smoother in this case. In case 11, the periodicity happens 
only for large flight distance \Sx\ with very small number 
of families. Note that the maximum of PDF is corre- 
sponding to the value \Sx\ of 2.2 which is just the length 
of a cell, and the PDF decays in power law. In this case it 
requires more collisions and takes more time for the par- 
ticles to escape a certain region. Thus the propagation 
is suppressed but is still of superdiffusion. The normal 
diffusion takes place when the particles are scattered by 
sufhciently large density of hyperbolic scatters (case IV) . 
Consequently, the strong chaos presents the trajectory 
of heat carriers with more aperiodicity. The PDF takes 
on its characteristic form which has a Gaussian tail as 
shown in the inset of Fig. |3Jd). 

Thus, the propagation modes are responsible for the 
diffusion behavior. The abundance of aperiodicity of tra- 
jectory is the characteristic of chaotic channel and may 
also play an crucial role in the normal diffusion. In other 
words, if the trajectory in a certain system emerges the 
aperiodicity due to some other mechanisms, such as in 
polygonal billiard gas model J^], the normal diffusion 
behaviors may happen. 



IV. THE CALCULATION OF TEMPERATURE 
FIELD 

We calculate the temperature field following the ap- 
proach proposed in Ref. The temperature of ith 
cell is defined by averaging the kinetic energy over all 
visits into the cell 



m 



(5) 



where tj denotes for the time spent within the cell in 
the jth visit, and m for the total number of visit. For 
sufhciently large m we expect a steady temperature pro- 
file, and this is indeed verified in our calculations for to- 
tally 10^° visits. The temperature profiles we obtained 
are plotted in Fig. |S1 It is worthwhile to point out 
that the steady temperature profiles between non-chaotic 
and chaotic system are quite different in thermodynam- 
ics limit, which is due to the different diffusion behaviors 
as shown in Fig. 0{c). As case I is non-chaotic and has 
uniform diffusion exponent the temperature profiles 
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keep almost the similar shape for different system sizes. 
At the two ends of channel there are large temperature 
jumps which play an important role in the Fourier trans- 
port and dynamics of the system |2 1| . These jumps arise 
from the boundary heat resistance which usually appears 
when there is a heat flux across the interface of the two 
adjacent materials. In case I and n which are chaotic 
and have asymptotically decreasing (3 versus time there 
also exists the boundary heat resistance. Unlike in case 
I, the temperature jump here is smaller and diminishes 
when the system size grows. For larger size is there al- 
most no temperature jump which corresponds to a nearly 
linear temperature profile. Both the larger system size 
A'' and the arc radius lead to the increase of chaos degree 
which is responsible for the decrease of diffusion expo- 
nent (3 1). In case IV, which is strong chaotic, the 
temperature profiles are almost linear for various system 
sizes we considered, corresponding to the normal heat 
conduction. 
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n-Lix) ^ 1 — X when /? = 1 (normal diffusion) and 
nL{x) ^ const when (3 — 2 (ballistic diffusion). The 
conservation of particle number requires 



Dl 



(6) 



where Dl are the diffusion coefhcient. Likewise, for a 
particle propagating from right to left, we have 



We assume Dl = and Dr — -^r 

2 /? ^ 1, the temperature is given by 



r; 



./3/2 



(7) 



Thus, if 



T{x) = 



TLnL{x) + TrHr^x) 



nL{x) + nFi{x) 
TlT^''^{1 - x)^ + t'I^'^Trx^ 



rf/V+T]^'"(l-a;)'' 



(8) 



where 7 = {2/ (3 - 1)(3^''^ . 




FIG. 6: (color on line) Numerical results of temperature pro- 
files for Tl = 1.0, Tr = 0.9 and sizes N = 20(solid), 40(dash), 
80(dot), 160(dash dot), 320(dash dot dot), 640(short dash) 
and 1280(short dot), respectively. The four panels refer to 
(1)7? = 0, /i = 1.0; (H) R = 0.001, h = 1.0; (H) R = 
0.848528, h = 1.0, and (IV) R = 0.848528, h = 0.27. The 
red lines correspond to the best fits for the numerical temper- 
ature profile at Ai' = 1280 with Eq. (|HJ, giving the analytical 
values 13 with 1.65, 1.34, 1.06 and 1.00 for above four cases 
respectively. 

We estimate the temperature profiles from the average 
point of view. Considering the incident particles from the 
left heat bath (where x — Q) propagating along the x-axis 
to the right end, we suppose that a reflecting boundary 
is placed at the origin of the x-axis and an absorbing one 
at the other end. When 2 ^ /? > 1, we assume that the 
mean density nL{x) of the particles at site x in the steady 
state is proportional to (1 — x)'* with 7 = {2/(3— 1)/?^^^, 
where we set x = i/N. Under this assumption, we have 



FIG. 7: (color on line) Numerical results of temperature 
profile in comparison to the analytical results. The nu- 
merical temperature profiles with R = 0.848528, A'' — 40; 
R = 0.4, iV = 40; i? = 0.2, = 40; i? = 0.1, AT = 80; 
R = 0.05, A' = 160 at /i = 1.0 and 7? = 0.848528, Af = 20 
ai h — 0.5 almost share the same shape. The red line is the 
plot of Eq.® with 13 = 1.25. 

As shown in Figs. |niand[7| the analytical results (in red 
lines) are in good agreement with the numerical ones for 
all the cases except those at the two ends of the channel 
for superdifusion cases. These deviations are likely due to 
the different boundary conditions we used. Furthermore, 
the value of (3, obtained by the best fits for the numeri- 
cal temperature profile at A^ = 1280 with Eq. ©, agrees 
with the simulating result greatly for aforementioned four 
cases. One can see clearly from Eq.® that temperature 
profiles are closely related to the diffusion exponent (3, 
namely, the case with smaller diffusion exponent tends 
to have smaller temperature jump. Accordingly, it is not 
unexpected that different chaotic cases may share the 
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same temperature profile if they have the identical dif- 
fusion exponent (3. As shown in Fig. [T] the case with 
smaller diffusion exponent requires smaller system size 
for achieving the same temperature profile. Moreover, 
our calculations show that the results of Eq. ((Sj) are con- 
sistent with the numerical ones even in larger tempera- 
ture gradient. Thus, the temperature profile is mostly 
dependent on the diffusion behavior which is remarkably 
affected by the finite-size effect for chaotic cases of su- 
perdiffusion. 

V. DISCUSSION AND CONCLUSION 

In summary, we have investigated the role that the 
chaos plays in the heat conduction by billiard gas chan- 
nel. We have demonstrated that the degree of dynamical 
chaos is enhanced by increasing the arc radius or the 
system size for chaotic channel, and the mass and heat 
transport behavior is significantly related to the degree 
of dynamical chaos of a channel. The stronger the chaos 
is, the closer to normal transport behaviors the model 
seems to be. Furthermore, our numerical results of two 
exponents a and [3 for both non-chaotic and chaotic cases 
when /3 ^ 1 satisfies the formula a = 2 — 2//3 [l^l- We 
also discussed the microscopic dynamics by the PDF of 
fiight distance in x-direction. It seems that aperiodicity 
of trajectory plays an important role in diffusion behav- 
ior. Finally, our results showed that the temperature 
jumps at both ends of the channel depend mostly on the 
diffusion property for both non-chaotic and chaotic chan- 
nels, and the finite-size effect is more crucial for chaotic 
ones. 



As is known that the billiard gas model is applicable 
for capturing the underlying dynamics of particles with- 
out interaction. It is therefore worthwhile to discuss the 
relation between our work and others in this field. 

Alonso et al. investigated ID Lorcntz gas model 
full of periodically distributed half circular scatters. By 
defining the heat conductivity and temperature field as 
statistical average over time on the hypothesis of local 
thermal equilibrium, the Fourier law holds in this case, 
and a linear gradient is given for quite small temperature 
difference . Our work starts from the same approach but 
different scatter geometry is taken into account. Thus 
it is not surprising that our work has some overlap with 
theirs in spirit. However, we pay much more attention 
to the role played by different degree of dynamical chaos 
in heat conduction. As a result, our intensive calcula- 
tions extended the results in Ref. and concluded 
that only sufficient strong chaos results in the normal 
diffusion, thus, the normal heat transport. 

Li et al. presented the dependence of heat con- 
ductivity on system size and the temperature profile in 
channel with zero Lyapunov exponent where the right tri- 
angle scatters are periodically distributed. In this case, 
the exponent stability leads to abnormal transport be- 
havior. Clearly, their result is the non-chaotic limit of 
our model. 
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